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In this Letter, we develop a microscopic theory to describe the close proximity between the insulating 
antiferromagnetic (AF) order and the d-wave superconducting (dSC) order in cuprates. We show 
that the cuprate ground states form a configuration of coherent pairing states consisting of extended 
singlet Cooper pairs and triplet tt pairs, which can simultaneously describe AF and dSC orders. 

PACS numbers: 74.20.-z, 74.20.Mn, 74. 25. Ha, 71.10.-w 



Discovery of high-T c superconductivity in copper- 
oxides ]l| reveals very attractive, but also extremely com- 
plicated, new phenomena in strongly correlated electron 
systems Q . One of the most striking phenomena is the 
close proximity between the Mott insulating antiferro- 
magnetic (AF) order and the d-wave superconducting 
(dSC) order. Theorists have tried various different mech- 
anisms, such as the theory of resonating valence bond 
(RVB) states of singlet pairs introduced by Anderson a 
decade ago || and the SO (5) unified theory of AF and 
dSC order parameters proposed by S. C. Zhang recently 
to explain this metal-insulator transition. However, a 
microscopic description of the AF to dSC transition has 
not yet been completed. 

It is believed that the underlying dynamics of strongly 
correlated electrons in cuprate superconductors can be 
modeled by the one-band 2D Hubband model or the t — J 
model HU . Unfortunately, the t — J model itself is diffi- 
cult to solve. Anderson's RVB state is a Mott insulating 
singlet at half-filling. However, it does not have AF or- 
dering. In Zhang's SO (5) theory, AF order parameters 
are a part of the building blocks. However, near half- 
filling, this AF order may not be Mott insulating, due 
to the lack of the constraint of no-doubly occupied sites 
which can destroy the SO (5) group structure. Gauge 
theories later developed in terms of slave bosons in the 
charge spin separation scenario of RVB states [||, and 
various microscopic tests of the SO (5) theory made in the 
past year J5J , certainly provided deeper understanding to 
high T c superconductivity, but the mechanism concern- 
ing the AF to dSC transition remains unsolved. 

Most likely, the low energy degrees of freedom in 
cuprates should still be dominated by electron pairs, al- 
though pairing structures may be different from the sin- 
glet Cooper pairs in BCS theory or the triplet pairs in He 3 
superfluid theory. Experiments confirmed the existence 
of the d x 2_ y 2 pairing symmetry in cuprate superconduc- 
tors. However, this does not imply that these electron 
pairs must all be spin singlet pairs, even though cuprate 
ground states are total spin singlet. Indeed, with d-wave 
singlet pairs alone, it is difficult to make a manifestation 
of the AF to dSC transition. 

In this Latter, we start with an insulating AF ordering 
of the undoped (half-filling) cuprate compound to ex- 



plore how superconductivity emerges in doped cuprates. 
We find that insulating AF and dSC orders can be si- 
multaneously described by a configuration of coherent 
pairing states consisting of extended singlet Cooper pairs 
and triplet 7r pairs. At half-filling, this configuration is 
dominated by the insulating AF ordering arisen from the 
mixing of singlet and triplet pairs. Upon doping with 
holes the d-wave pairing gap appears, and above a cer- 
tain level cuprates become superconducting. In this de- 
scription, the nature of no-doubly occupied sites and the 
present of triplet 7r pairs play an important role. 

We shall begin with the formation of electron pairs in 
strongly correlated systems. The square lattice structure 
of layered copper-oxides allows us to restrict the problem 
in the reduced (half) first Brillouin zone in momentum 
space. Then the cuprate ground states may be obtained 
by projecting out all states of no-doubly occupied sites 
from the following generalized pairing state $}: 

= Il k cx P {^( k ) c h c -U + »»( k )4+QT c -k+Qi 
+ %(k)4 T c! k+Qi + ?74(k)cj ci cL k+QT 
+ %(k)c kT c! k+QT + %(k)c ki c T _ k+Qi 
-H.c}|0). (1) 

where the complex parameters ?7i(k) are generally link- 
dependent pairing wave functions. The parity symmetry 
in the x — y plane requires ^(k) = k). In the gen- 
eralized coherent state theory ||, $} is a product (over 
k) of the local SO(8)/U(4) coherent pairing states Q. 

In fact, Eq. ([j]) is also a multi-pair generalization of 
BCS pairing states or RVB states (after projection). Ex- 
plicitly, we may rewrite ?7i(k) = rj s (\t) + i]d(k), i]2(k) = 
r] s (k) - r/ d (k), r? 3 (k) = ?/V (k) + ^(k), ry 4 (k) = r/V (k) - 
77,, (k), and i] 5 (k) = r) n+ (k), r/ 6 (k) = i] 7r _(k). Then, 
Eq. (Q) consists of all electron pairs concerned in the 
study of superconductivity. These are the ordinary s- 
wave Cooper pairs given by the ?7s(k) term, the ex- 
tended (i.e. valence bond) singlet Cooper pairs by the 
77d(k) terms [including the extended s-wave, d-wave and 
s + id pairs in terms of ry ( j(k) = ^^(k), r/ k d(k) and 
r/ k +id (cos k x + i cos k y ) respectively], the quasi-spin rj 
pairs by the f/V^k) term which generate quasispin SU(2) 
transformations, the singlet p-wave pairs [if one lets 
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7777 (k) = ^(sinfca; ± sinfcj,)], and finally the triplet it 
pairs by JyV^k) = ry^c?(k), here 7(k) = cosk x + cosfc y 
and d(k) = cos k x — cos k y . 

At half filling, without any loss of generality, we set 
r] s (k) = r] v (k) = 0, i.e., ?7i(k) = -77 2 (k) and %(k) = 
774 (k) for the AF ordered phase. In addition, the cuprate 
ground states must also be spin singlet due to the spin 
rotational symmetry which implies ($|S ,C *|$) = 0. This 
can be satisfied if 775 (k) = — ^(k). These restrictions 
on pairing wave functions can be justified by explicitly 
calculating the AF order parameters from (|l|): 

m z = ~2jy^2 k [( z lk - ^2k)(>3k + z 4k) + H.C , 

m + = m_ = ^Xl k [( Zlk ~~ Z2k )( z 5k - z k) + H - c > ( 2 ) 

where m z and m± are the average staggered magneti- 
zations in the z direction and the x — y plane, respec- 
tively, N is the total number of lattice sites, and the 
parameters Zjk are elements of the 4x4 block matrix 
in the canonical Bogoliubov transformations induced by 
(0): = sin \J rfv/Vv^V with 



7?(k) 



(k) 775 (k) 773 (k) 

-77! (k) 774 (k) 776 (k) 

-775 (k) -774 (k) 772 (k) 

-773 (k) -776 (k) -772 (k) 



(3) 



The undoped ground state maximizes the staggered mag- 
netization which corresponds to 



Zlk 



-Z 2 k — Zdk 7 z 3k — ^4k , 2 5 k 



(4) 



This is the same as letting 771 (k) = — 772 (k), 773 (k) = 774 (k) 
and 775 (k) = — 776 (k). These restrictions on the pairing 
wave functions are also necessary for the manifestation of 
the extended (including <i-wave) pairing symmetry of va- 
lence bonds, because Eq. (Q) implies 77* (k + Q) = — ?7i(k). 

In the above calculation, we have not imposed the con- 
straint of no-doubly occupied sites. However, the result 
can be further varified under this constraint. Here, in- 
stead of using Gutzwiller projector Pq to remove doubly 
occupied sites from |$), it is equivalent to require that 1$) 
must satisfy the constraint ($| X)i( n jT n «l)l < ^ > ) = 0> wnere 
summation to i is over the lattice sites. This is because 
the operator n^riii is positive definite. In the mean- field 
approximation of (|lj), we obtain the global constraint: 



^-^ + |A s | 2 + |A,| 2 = 0, 
where n = ($|t7 j |<I>} is the total electron number: 



(5) 



(0) 



and 5 the fractional doping of holes, m s denotes the mag- 
nitude of the long-range AF order parameter determined 



from (||) by to 2 = to 2 + m + m-, and A s and are av- 
eraged order parameters of the ordinary s-wave Cooper 
pairs and the quasispin 77 pairs which vanish under (Q): 



A. 



0. Then, Eq. (§) reduces to 



4AT2- 



(7) 



0.5 which is the 



At half- filling (5 = 0), it gives m s 
same result as in Neel states. 

Furthermore, because of the spin rotational symmetry, 
without loss of generality we can define 2:3k = z^k cos 2#k 
and 2:5k = z,rk sin 2#k- Thus, from (^), we find that #k = 
9 for all k pi and 



Zirk = z d k or Z Tk = -Zdk- 



(8) 



This indicates that under the constraint of no-doubly oc- 
cupied sites, the extended singlet Cooper pairs and triplet 
7r pairs in |<I>) have equal occupied probabilities. Also, 
Eqs. (||) and (||) reduces to 
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= 8^'|z dk | 2 = iV(l-,J) (9) 



and m z = m s cos2#, m + = to_ = m s sin 29. Here m s 
and n can be written only in terms of the singlet pairing 
wave function z^k because of the constraint of no-doubly 
occupied sites (||). At half- filling, electrons are uniformly 
distributed. Then Eq. (H) gives 



\Zdk\ 



1/2. 



(10) 



This is indeed the local constraint of no-doubly occupied 
sites at half-filling. One can show that ($|(fi — n) 2 |(f>) = 



under (g) and (|10|). 

The above rigorous results are derived only from the 
kinematics of the coherent pairing state |$). Dynam- 
ically, we may determine the ground state from the 
t — J model. The i-term is a hoping Hamiltonian H t = 
Ek,o- £ ( k ) c L c ^' where e(k) = -2t(cosk x +cosk y ). Its 
expectation value is 



(Ht) 



2^ k e(k)(|z lk | 2 



z 2k 



(11) 



which vanishes by (^), as expected from the constraint 
of no-doubly occupied sites at half-filling. Under the re- 
striction (jj) , the expectation value of the J-term [Hj = 

J J2( t ,j)( S i ■ S j - \ n i n j)] is s iven b y 



(Hj) 



-(|A d | 



A d | 2 + |A C 



2m 2 



X 



JN V 1- " 1 ' ' — " s ' 2A 2 /' 

where the averaged <i-wave pairing order parameter 



(12) 



z£w± + z v w" 



N 

^^(k)z dk v / l-4|z d k| 2 



(13) 
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and the extended s-wave pairing order parameter A es 
has the same form as ( |l3| ) by replacing d(k) with 7(k). 

The parameters z^ = z^w ± z^k and = yjl — \z^\ 2 . 
Note that although it can be written in terms of z d k in the 
second equality of ( |l3| ) [because of the constraint (g)], 
receives contributions from both the singlet and triplet 
pairs. As a result, A^ behaves very different from that 
in pure singlet pair states. At half-filling, (Jl^) ensures 
Ad = A es = 0. Then (Hj) = —JN which again gives 
the same result as in Neel states. 

These results are very remarkable. It shows that un- 
der constraints (|J) and (||), the coherent pairing state 
|$) is insulating AF ordering at half-filling. It consists 
of extended singlet Cooper pairs and triplet tt pairs with 
equal pair numbers. The singlet pairs are indeed Ander- 
son's RVB pairs ||, while the triplet tt pairs are the tt 
operators of Demler and Zhang Both the AF and 

dSC order parameters are functions of extended pairing 
wave functions (including the extended s-wave, d-wave 
and s + id pairing functions Jl2|) of singlets and triplets. 
However, because of the constraint of no-doubly occupied 
sites, the pairing order parameter vanishes while the AF 
order parameter is the same as in Neel states. 

Meanwhile, although the configuration of coherent 
pairing states contains the same AF properties as the 
Neel configuration, these two configurations are differ- 
ent. Electrons in Neel states are not paired and Neel 
configuration itself is only valid for half-filling. However, 
|$) is defined for any doping and it is built with extended 
(valence bond) pairs. Thus, |$) provides a natural path 
to explore how cuprates undergo the transition from an 
insulating AF order to a dSC order after dopings. 

Upon doping with holes, it is commonly believed that 
hopings rapidly distroy the AF ordering. Eq. (|ll|) shows 
that the lowest energy gained from the t-hoping corre- 
sponds to z 2 k = 0, which reduces AF order parameter by 
a half, as shown by (||). Apparently this is a nice result 
as expected from hoping dynamics. However, z 2 k = 
leads to a non- vanishing A s . A non- vanishing A s en- 
hances the AF order parameter m s by the constraint (||) , 
which is contrary to (||) for Z2k = 0. Indeed, the ordi- 
nary Cooper pairs are on-site pairs which are excluded 
by the constraint of no-doubly occupied sites. In terms of 
energy scales, the on-site repulsive Coulomb interaction 
is much stronger than hopings, a non-vanishing A s con- 
tributes a large positive Coulomb energy to the ground 
states, which is certainly unfavored. Indeed, any net con- 
tribution from the i-hoping in mean-field theory requires 
\ z ik\ 2 l z 2k| 2 7 which will break the extended singlet 
pairing symmetry. Hence, doped ground states must still 
be restricted by (||). 

As a result, after dopings, |$) still consists of the ex- 
tended singlet Cooper pairs and triplet tt pairs. All equa- 
tions derived above are still valid except for @. Eq. (|) 
shows that the AF order parameter m s is decreased lin- 
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FIG. 1. Gap order parameters via dopings at zero temper- 
ature. The (top) solid line is for t'/J = 0, and the middle and 
bottom lines are for t'/J — 0.2 and 0.3, respectively. 

early with doping S. By further including the quantum 
fluctuation of spin- wave excitations, m s is suppressed 
from 0.5 to 0.3 at half-filling p3[ , and will be suppressed 
as well upon dopings. While, the pairing order parame- 
ter Aw [given by Eq. @] is zero at half-filling. Then it 
is increased after dopings so that superconducting states 
will emerge at a certain level of dopings. 

To be more explicit, we may take again the t — J model 
as an example. When all electrons in ground states 
are paired, the t-term vanishes by the pairing symme- 
try £ik = — -Z2k for any doing. However, the next order 
hoping (t'-term) has a non-zero expectation value in |$). 
The ground states can be determined by minimizing the 
t 1 — J Hamiltonian, which leads to z^k — ► and 
the gap equation at zero-temperature, 

A k = — yVkk- — . (14) 

where A k = A d d(k) + A es7 (k), U k k< = J[d(k)d(k') + 
7 (k) 7 (k')] and £ k = { J 2 Al + [e(k) - fx -2(1 - 8) J] 2 } 1 ' 2 
with e(k) = — 4i' cos fez cos fcy. The fixed electron num- 
ber gives -|Ek £(k) ^ (1 ^^ )J = 1-2(5, and fi is the 
chemical potential. 

The solutions of (|lj) show that the d-wave gap order 
parameter appears after dopings for S < 0.5 with max- 
imum peak A d ~ 0.07 ~ 0.10 at S ~ 0.15 - 0.20 (the 
typical optimal doping region) for t'/J = 0.30 ~ 0.20. 
This is in good agreement with the experimental obser- 
vations of d-wave superconducting states in cuprates. It 
is also striking that the extended s-wave superconducting 
states only emerge in overdoped region of S > 0.5. The 
separation of the d-wave states in optimal dopings from 
the extended s-wave states in overdopings is controlled 
by the £'-term. Here t' must be positive. For a negative 
t' , the ordering of A^ and A es in terms of 5 will be ex- 
changed. If we let t' = 0, then A^ = A es which have the 
maximum value at doping S = 0.5. This corresponds to 
the symmetry limit of Zhang's SO (5) theory [see ([l2h], 
although generally the pairing wave functions determined 
here do not form a rigorous global SO (5) group structure 
fuj . The numerical results are plotted in Fig. 1. 

The above results are very different from the d-wave 
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singlet pair theory. The d-wave singlet pair states corre- 
spond to the special case of (|l]) with 771 (k) = — 772 (k) = 
i]d(k) and other r]i(k) = (i =3 to 6). Then from (0), 
(0) and (|l|), we have 

m s = , n = 4V ' \z dk \ 2 = 7V(1 - J) 

* — 'k 

A d = ^V' d(k)z dkv /l - | 2dk | 2 . (15) 
iV ^ — 'k 

Hence, with singlet pairs alone, no AF order exists. The 
order parameter has a maximun value near half-filling 
(zxk — \/l/2), then it is slowly decreased with dopings 
fjjfl . Of course, here the constraint of no-doubly occu- 
pied sites cannot been imposed analytically by (hj) so that 
auxiliary-boson variables must be introduced |q] . 

Formally, (|l5| ) and (|l3|) show that the magnitude of 
in pure singlet pair states is about twice of that in mixed 
states of singlet and triplet pairs. In other words, pure 
singlet pair states have lower gap energies than mixed 
pair states. However, mixed pair states maximize the 
AF energy which cannot be taken into account in pure 
singlet pair states [see Eq. ( |l2] ) in the t — J model, for an 
example] . Because of the AF energy, mixed pair states 
have indeed a lower ground state energy than pure singlet 
pair states in the mean-field theory. 

Summarizing the above discussions, we show that un- 
derlying pairing states undergoing an insulting AF to 
dSC transition can be determined by the coherent state 
$) of extended (valence bond) singlet Cooper pairs 
mixed with triplet n pairs. These coherent pairing states 
are unconventional, and can be regarded as a nontriv- 
ial generalization of Anderson's singlet pair RVB states. 
Advantage of the coherent pairing state $) is that it con- 
tains an intrinsic AF order which is Mott insulating at 
and near half-filling, while in optimal dopings, it has the 
d-wave SC order. Therefore, it can microscopically unify 
the insulating AF and dSC orders in the same configura- 
tion. This is for the first time to analytically show that a 
pairing state can simultaneously describe both the AF or- 
der and the dSC order. Dynamically, we show that due to 
the extended singlet pairing symmetry, the contribution 
from the t-hoping is, surprisingly, cancelled in (classical) 
ground states. The exchange interaction ( J-term) mani- 
fests a symmetry limit of Zhang's SO (5) theory. While, 
the insulating AF to dSC transition under dopings that 
is comparable with experimental observations can be de- 
termined by the next nearest-neighbor hopings (with a 
positive but weak t' hoping). 

Furthermore, the canonical transformation structure 
of (^) generates a quasiparticle picture, which allows us 
to easily determine normal state properties of cuprates 
as well as the phase diagram of high T c superconduc- 
tivity in the T — S plane [jf5j. Also note that cuprate 
ground states determined here by minimization are classi- 
cal states. In these classical ground states, z^ — > |^dk|e^ 
and 6*k — » 9 so that the configuration manifests the global 



SU spin (2)xU gaugc (l) degeneracy. This SU(2)xU(l) sym- 
metry will be spontaneously broken at quantum level. 
The quantum fluctuations of Zdk. and #k describe pairing 
excitations and spin- wave excitations (Goldstone modes) 
which further lower energies of the ground states over 
all dopings. Therefore, these fluctuations play a very 
important role in the determination of dynamical prop- 
erties of electrons in cuprates. A typical evidence is the 
AF magnetic ordering which is too strong in classical 
ground states [see (0)]. However, it can be largely sup- 
pressed (and should be cancelled in optimal dopings) by 
the quantum fluctuation of spin- wave excitations |Hj| . In 
separate publications jL6), based on the generalized co- 
herent state theory || we shall use |$) to develop a theory 
of quantum fluctuations to detail these features, and then 
to describe dynamical and thermal properties of high T c 
superconductivity in cuprates. 
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